Bose-Hubbard Models Coupled to Cavity Light Fields 
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Recent experiments on strongly coupled cavity quantum electrodynamics present new directions in 
"matter-light" systems. Following on from our previous work [Phys. Rev. Lett. 102, 135301 (2009)] 
we investigate Bose-Hubbard models coupled to a cavity light field. We discuss the emergence of 
photo-excitations or "polaritons" within the Mott phase, and obtain the complete variational phase 
diagram. Exploiting connections to the superradiance transition in the Dicke model we discuss 
the nature of polariton condensation within this novel state. Incorporating the effects of carrier 
superfluidity, we identify a first order transition between the superradiant Mott phase and the single 
component atomic superfluid. The overall predictions of mean field theory are in excellent agreement 
with exact diagonalization and we provide details of superfluid fractions, density fluctuations, and 
finite size effects. We highlight connections to recent work on coupled cavity arrays. 

PACS numbers: 03.75. Mn, 03.75. Hh, 67.85.-d, 05.30.Jp 



I. INTRODUCTION 

Over the last few years there has been tremendous ac- 
tivity in studying the coherent interaction of matter and 
radiation in a rich variety of "matter-light" systems. Re- 
cent experiments have combined cavity quantum electro- 
dynamics (cavity QED) with cold atomic gases, and allow 
access to the strongly coupled regime. This has led to pio- 
neering work on Bosc-Einstcin condensates (BECs) in ul- 
tra high finesse optical cavities [l[ , and with optical fibers 
on atom chips Q • It has also stimulated advances in cav- 
ity opto- mechanics and condensate dynamics More 
recently, strong matter-light coupling has been achieved 
for ion crystals with potential applications in quan- 

tum information processing. These developments offer a 
wealth of possibilities, at the interface between quantum 
optics, cold atoms and condensed matter physics. The 
light field serves not only as a probe of the many-body 
system Q, but may also support interesting cavity me- 
diated phenomena and phases. It may further provide 
routes to simulate strongly correlated quantum systems, 
with proposals based on coupled cavity arrays [8142 31] and 
nonlinear optical fibers [24| . 

Allied advances in solid state devices include cavity 
QED experiments with superconducting qubits in mi- 
crowave resonators [25j . This has provided clean real- 
izations of the paradigmatic Jaynes-Cummings f2fj| and 
Dicke models [27|, describing two-level systems coupled 
to radiation. It has also led to remarkable observations of 
the Lamb shift [28| . This is complemented by the quest 
for pol ariton condensates in semiconductor microcavities 
[29l434| , where the hybridization of an exciton and a pho- 
ton yields low effective mass polaritons. This offers the 
prospect of higher transition temperatures than for ex- 
citon BEC, and gives access to coherence properties via 
the cavity light field. 

Motivated by this broad spectrum of activity, we ex- 
amine the impact of cavity radiation on bosonic Hubbard 
models 35 1. Following on from our previous work [36j . 



we focus on a two-band model in which photons induce 
transitions between two internal states or Bloch bands. 
This is a natural generalization of the much studied two- 
level systems coupled to radiation, and may serve as a 
useful paradigm in other contexts. In Ref. [3fj| we dis- 
cussed the interplay of Mott physics, photo-excitation, 
and Bosc condensation promoted by carrier itinerancy. 
In particular, we provided evidence for a novel Mott 
phase with photo-excitations analogous to polaritons. In 
this work we study this problem in more detail, with em- 
phasis on the nature of the polariton condensate. We 
also highlight connections to coupled cavity arrays de- 
scribed by the Jaynes-Cummings-Hubbard model and 
its variants 

HMMSl- Additional directions in cold 
atoms include recent work on excitons (37j . generalized 
Dicke models K58( , and light propagation in atomic Mott 
insulators |39l. |40|| . 

The outline of this paper is as follows. We begin in 
section |TT] with an introduction to the two-band Bose- 
Hubbard model coupled to quantum light [3(| . In section 
IIIII we discuss the zero hopping limit of this model and 
anchor the phase diagram. We study the evolution with 
the strength of the matter-light coupling and highlight 
the connections to the Dicke model and the superradi- 
ance transition [4ll - |44| . In section IIVI we use a varia- 
tional approach to obtain the overall phase diagram in 
the presence of hopping. We corroborate our findings in 
section [V] using numerical simulations. We discuss the 
relation to the bosonic BEC-BCS crossover and to other 
problems of current interest in section I VI I We conclude 
in section IVTT1 and provide technical appendices. 



II. MODEL 

Let us consider a two-band Bose-Hubbard model cou- 
pled to the quantum light field of an optical cavity in the 
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FIG. 1: (Color online). Heuristic diagram showing photo- 
excitations above the filled bosonic Mott state. In analogy 
with fermionic band insulators, the coherent superposition of 
a particle-hole pair and a photon will be referred to as a "po- 
lariton". Within an equilibrium framework these polaritons 
may Bose condense. 



rotating wave approximation [36[ : 
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where i labels lattice sites, and a — a, b are bosons 
obeying canonical commutation relations [a, , aj] = <5,j. 
These might be states of different orbital or spin angular 
momentum, and we assume the cavity radiation field, if), 
may induce photo-excitations between these states; see 
Fig. [TJ Here, e a effects the band splitting ujq = e& — e , 
and U aa i are interactions, where :: indicates normal or- 
dering. This yields : ni a rii a := rii a (ni a — I) for like 
species and : ni a rii a i := ni a rii a > for distinct species. J a 
are nearest neighbor hopping parameters, and lu is the 
frequency of the cavity mode. For simplicity we con- 
sider just a single mode which couples uniformly to the 
bands. The coupling g is the strength of the matter- 
light interaction. In view of the box normalization of the 
photon field, the dipole coupling strength is proportional 
to 1/VK, where V is the volume of the cavity. With a 
fixed density of lattice sites, p = N/V, it is convenient to 
denote g = g/VN, where N is the total number of lat- 
tice sites. We work in units where the half-band splitting 
wo/2 = (e 5 - e„)/2 = 1. 

An important feature of the Hamiltonian (pj is that the 
individual atom and photon numbers are not conserved 
due to the matter-light interaction. However, the total 
number of atomic carriers 



(2) 



and the total number of photo-excitations 

N 2 = ^+^( n i- n i +1 ^ ( 3 ) 

i 

are conserved and commute with Hq. The latter counts 
the total number of photons plus particle-hole pairs, and 
we refer to these composite excitations as "polaritons" 



— see Fig. [TJ These conservation laws reflect the global 
U(l) x U(l) symmetry of H such that 



a -> e"V, b -> e llp b, ip -> e 



-i(ii-ip) 



(4) 



where i),^eR. In general this symmetry involves mix- 
ing between the matter-light sectors. Moreover, it also 
allows for the simultaneous coexistence of Mott behav- 
ior and condensation, corresponding to an unbroken U(l) 
and broken U(l) symmetry respectively. This symmetry 
will have a direct manifestation in the phase diagram, 
and will suggest implications for other multicomponcnt 
problems. We will work in the grand canonical ensemble: 



H = H - niNi - n 2 N 2 . 



(5) 



We begin by assuming that a are strongly interacting 
hardcore bosons, U aa — > 00, and that 6 are sufficiently 
dilute so that we may neglect their interactions, Ubb = 
0. We will also start with U a t, = 0, before discussing 
departures from these conditions. 



III. ZERO HOPPING LIMIT 

Before embarking on a detailed examination of the 
model (UJ, we investigate the zero hop ping li mit. As 
in the single-band Bose-Hubbard model [43, Sa] this an- 
chors the topology of the general phase diagram. In the 
present case this is particularly informative since the zero 
hopping phase diagram evolves with the matter-light 
coupling, g. In addition, the global photon mode cou- 
ples all the sites, even in this zero hopping limit. To gain 
a handle on this reduced Hamiltonian, we find it conve- 
nient to develop two complementary approaches. These 
help illuminate different aspects of the more general itin- 
erant problem, and provide a platform for extensions. In 
scction llll Al we begin with a variational approach using a 
coherent state ansatz for the photons. This will enable us 
to derive an effective single site Jaynes-Cummings model 
471 ] for the a, b bosons, and proceed with the minimum of 
technical input. In section fill Bl we instead map the a, b 
bosons on to effective spins. This yields the paradigmatic 
Dicke model [4l| , describing many spins coupled to radia- 
tion. The Jaynes-Cummings and Dicke models are famil- 
iar in atomic physics and quantum optics and we provide 
a brief overview of these closely related Hamiltonians in 
Appendix [A] Both approaches yield equivalent results 
and indicate a novel quantum phase transition occurring 
within the lowest Mott lobe [36| . This takes place from a 
conventional Mott state with no photons, to a state with 
a non-vanishing population of photo-excitations, as the 
strength of the matter-light coupling is increased. This 
quantum phase transition coincides with the well known 
superradiance transition in the Dicke model (4lH44l | , and 
we interpret the results in this framework. The combina- 
tion of perspectives will be useful in the generalization to 
the itinerant problem discussed in section ITVl For other 
proposals of Dicke models and sup erradiance transitions 
in cold atomic gases see Refs. j48H50| . 
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A. Variational Approach 

In the zero hopping and U aa ~> oo limit, the Hamilto- 
nian ([5]) may be written in the form 



i i 



(6) 



where the site occupancy of a-atoms is limited to 0,1. 
We absorb the chemical potentials in to the definitions 



e a = e a - Mi + M2/2, 
e 6 = e 6 - Mi - M2/2, 

LU = UJ — /i2- 



(7) 



This bosonic Hamiltonian is analogous to the fermionic 
problem of localized excitons in a microcavity [Fll . l52j . 
An instructive way to analyze the problem is to consider 
a coherent state for the cavity light field: 



(8) 



where 7 is a variational parameter to be determined. As 
usual, this is an eigenstate of the annihilation operator, 
ip, with eigenvalue 7. The coherent state assumption is 
expected to become exact in the thermodynamic limit, 
and this is borne out by the complementary approach in 
section MI Bl The expectation value 



(9) 



yields an effective single site problem for the a and b 
atoms which is readily diagonalized: 



+ e b n b +W7 2 + gjfta + alb), (10) 



where Cu = to/N. This is a significant merit of the current 
approach, since it is directly tractable without technical 
input. In view of the hardcore constraint on the a atoms, 
this describes a single two-level system coupled to an ef- 
fective "radiation field" of b atoms. This paradigm is 
described by the much studied Jaynes-Cummings model 
[UIH^], discussed in Appendix lAl The eigenstates of (fT0|) 
are superpositions in the upper and lower bands. Focus- 
ing on Mott states with total occupation n a + n b = n, we 
consider admixtures of |0,n) and |l,n— 1) in the \n a ,n b ) 
basis. The lowest eigenstate has energy 



E- 



W7 



ne b 



wo/2 - v /^/4 



+ g 2 7 2 n 



111) 



where cjo = Q> — ?a is the effective band splitting. The 
nonlinear dependence on ^Jn is a notable feature of the 
Jaynes-Cummings eigenstates and has recently been seen 
in circuit QED experiments [2(| . (Analogous dependence 
is also seen in BECs [HQ as a function of the number of 



atoms in the cavity.) Minimizing with respect to 7, one 
obtains the variational self-consistency condition 
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2 7 



g 2 n 



7 2 7 2 n 



0. 



(12) 



Depending on the parameters, one may therefore obtain 
cither the trivial solution 7 = 0, corresponding to zero 
photon occupancy, or the non-trivial solution 



7 2 

/var 



1 



2 

9 n 



g 2 n 



(13) 



corresponding to a finite photon occupancy. This latter 
solution is supported in the region where 7 2 ar > 0, or 
when the matter-light coupling exceeds the critical value 

Qc = 9c\fn = \/^o, (14) 

where for simplicity we consider ujq > 0. As we shall 
discuss in section IIIIB1 this onset of the photon field 
corresponds to the well known supcrradiancc transition 
in the Dicke model jriTli^ l . Indeed, it is readily seen from 
equation (| 13[) . that for Q > Q c , the photon field has the 
characteristic variation 



N 



1 

4^ 



Qt 



(15) 



where Q = g~\fn. See for example Table 1 of Ref. [54 1. 
We shall give an alternative derivation of these results 
in section IIII B[ This complementary approach becomes 
asymptotically exact in the thermodynamic limit and 
helps justify the coherent state ansatz ([5]). 

Knowledge of the eigenvalues (fTTj) . together with the 
variational consistency condition (|12l) . enables one to 
construct the zero hopping phase diagram depicted in 
Fig. [5] The lower boundary consists of two segments. 
The first corresponds to the transition from the vacuum 
to the Mott state with no photons, and is given by the 
condition E± (7 = 0) < or /ii > e a + M2/2. The second 
segment corresponds to the transition in to the photon 
rich Mott state and occurs when £^(7 var ) < where 



j?-, \ ~ ^0 n (g A + gt 
E n {l™)=ne b ----{- ir - 



Explicitly this locus is given by 



in 



> 



a, 



1 

4tD 



+ U) 2 {UJQ — /Lt 2 ) 5 



(16) 



(17) 



Noting that g c depends on n according to ([T4|) we see that 
S,7(7var) -> n(e b -g 2 /4u>) as n -> 00. For e b -g 2 /4ui < 0, 
or chemical potentials satisfying 



Ml > e b 



/'2 



4(cj - M2) ' 



(18) 



it is energetically favorable to macroscopically populate 
the b levels. This corresponds to the upper boundary in 
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FIG. 2: (Color online). Zero hopping phase diagram in the 
U aa — > oo limit with Utb = U a t = 0. We set tb = — e a = g = 1 
and w = 1 corresponding to u < ojq, where u>o = £b — e Q . 
The vertical line, g = g c , is the superradiance transition in 
the Dicke model, and separates a Mott insulator (MI) with 
n a + rib = 1 and (ip^4>) = 0, from a superradiant Mott in- 
sulator (SRMI) with (tp^ip) 0. Outside of these regions 
are the vacuum state, and the unstable region corresponding 
to macroscopic population of the b states. Whilst the total 
density is fixed within both Mott phases, n a + rib = 1, the 
individual a and b populations vary in the superradiant phase 
as shown in Fig. [3] Inset: eb = —e a = g = 1 and u) = 3. 
For ui > uio the upper and lower boundaries may cross and 
terminate the lobe, since for /12 > ^0 the lowest stable state 
is the vacuum; see Appendix [D] for a derivation of this. 




0.25 0.5 0.75 



FIG. 3: (Color online). Onset of magnetization A4 = 
(rib — n a )/2 at the superradiance transition, traversing a 
line through Fig. [2] with /ii = 0. Whilst the total density 
n a + rib = 1 remains pinned throughout the entire Mott re- 
gion, photo-excitation promotes atoms between the bands in 
the superradiant phase. 



Gc- In Fig. |3]we plot the continuous onset of population 
imbalance described by equation (|20|) . In section MI Bl 
we will see how these results emerge from the exactly 
solvable Dicke model. 



B. Dicke Model 



Fig. [2] Indeed, in the absence of the matter-light cou- 
pling equation (fT5)l yields e& < and is naturally associ- 
ated with population of the b states; we provide the cor- 
responding .g = phase diagram for comparison in Ap- 
pendix [B] In contrast to the classical light case (where i/j 
is treated as a fixed c-number) coupling to the fluctuating 
quantum field, ip, eliminates the higher Mott lobes cor- 
responding to integer increases in the ^-populations; the 
higher lobes corresponding to the a-particles have been 
explicitly eliminated by the hardcore constraint. We con- 
firm this reduction analytically in Appendix [Dl 

Accompanying the superradiance transition is a change 
in the relative atomic populations or "magnetization" 



M = - ((n b ) - (n » , 



(19) 



corresponding to photo-excitation into the upper band. 
Indeed, the possibility of such magnetic phase transitions 
is a strong motivation for studying multicomponent prob- 
lems, even without matter-light coupling [5TjI [56l|. The 
magnetization (fTT))) is readily computed from the Jaynes- 
Cummings cigenstates (|A6|) with n = 1. It also serves as 
an order parameter for this continuous transition: 



M = 



Q < Gc, 



-4 % ; Q>Qc 



(20) 



This corresponds to (n a ) = 1 and (n b ) 
Gc, and a non-trivial imbalance, (rib) — 



for g 
for G 



As in the fermionic cases considered elsewhere [5l|, [f52j , 
an alternative way to view the zero hopping Hamiltonian 
© is as an effective spin-boson model. Within the sub- 
space of fixed density, n a +n b = n, we introduce effective 
spins for a priori possible Mott lobes 

|^) = |l,n-l), |ft) = |0,n>, (21) 

where we denote bosonic states as \n a , rib) ■ The operators 

S+= b -^, S~ = ^, S* = \-n a , (22) 

form a representation of su(2) on this restricted Hilbert 
space. In the lowest Mott lobe with n = 1 this reduces 
to the usual Schwinger boson construction [57j |. In the 
representation (|2"2")) the Hamiltonian © becomes 



JV 



N 



H = Q ^Sf 



where G = and c n = N(nib — Qg /2) . This is the 

much studied spin-1/2 Dicke model [41(, describing N 
two-level systems coupled to radiation; see Appendi x [XI 
for a brief review. This model is integrable [U H, |58| | 
and in the thermodynamic limit it has a quantum phase 
transition to a so-called superradiant phase when |43 - [44| 



Gc = Gc^N = y/, 



UJUJo, 



(24) 
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in agreement with our previous result (fT4|) . In this con- 
text, the term superradiance indicates the onset of many- 
body or cooperative effects involving the photon coupled 
to many atoms. This mapping not only helps justify our 
variational approach, but will also provide insights into 
polariton condensation and matter-light coherence at the 
superradiance transition. 

The thermodynamic limit of the Dicke model ([23]) may 
be analyzed using collective spin operators [59j 



(25) 



where we exploit site independence of the global photon 
field, and N = 2S plays the role of a large spin. This mo- 
tivates an asymptotically exact semiclassical treatment 
based on the Holstein-Primakoff transformation [59146 1| 



^(25- 



UV2 



c' c 



J" 

J- = (25-c t c) 1 / 2 c, 
J z = Jc - S, 



(26) 



where c is a canonical boson. The superradiance tran- 
sition is associated with condensation of this auxiliary 
boson. Most crucially, this is related to condensation of 
polaritons, and not the a, b bosons themselves. The Dicke 
model (12"31 becomes 



H = u>o (c f c - S) + atyty 

+ _J= ^t (25 _ C t c) i/2 C + c t(25 _ c t c )i/2^ 



where we drop the 'constant' c n . Within a semiclassi- 
cal 1/S expansion around the thermodynamic limit [59[ 
we may proceed by introducing coherent states for the 
photon and auxilliary boson 



|o>, 10 = 



|0), (28) 



where 7 and £ are classical c-numbcrs. The energy den- 
sity reads 

£ = (H)/N = £ (C 2 -l/2)+£7 2 + 2S7C(l-C 2 ) 1/2 , (29) 
where 7 = j/y/N, C = C/V^V- Minimizing over 7 gives 

_ eca-c 2 ) 1/2 



Substituting back in (f2"§[) yields 



U \ U! J 2 



(30) 



(31) 



and we acquire an expectation value (c) 7^ when the 
quadratic term becomes negative. This corresponds to 
the superradiance quantum phase transition at Q c = 



\ZCjCjq in agreement with the previous results. Minimiz- 
ing ([31]) over ( 2 and substituting into ([30]) one obtains 



C 2 



2Q 2 



7 



AQ 2 G 2 



Q>Qc 



(32) 



The onset of the photon field agrees with equation (|15p . 
In addition, the expectation value ( J + ) = (,(2S — Q 2 ) 1 ' 2 
tracks the condensation of the Holstein-Primakoff boson, 
c, via its coherent state parameter Q: 



(J- 



N 



c(i-c 2 ) 1/2 = -^. 



W7 



(33) 



These results show that the superradiance transition is 
accompanied by condensation of the Holstein-Primakoff 
boson, (c) 7^ 0. As we will discuss in section IIVB1 
in the two-band Bose-Hubbard model coupled to light, 
this corresponds to condensation of the polaritonic bilin- 
ear (b^a) above the Mott background, and not the a, b 
bosons themselves. In addition to the onset of photons 
described by equations (|15|) and (|32[) . the magnetization 
A4 = (J z )/N reproduces the previous results. 

The correspondence between the variational approach 
outlined in section IIII Al and the Dicke model analysis 
is clearly encouraging. In section ITvl we shall extend the 
variational approach to include the important effects of 
itinerancy and carrier superfluidity. 



IV. VARIATIONAL APPROACH FOR 
HARDCORE ATOMS 

A. Phase Diagram 

Having confirmed a zero hopping Mott phase, with 
n a + n>b = 1, we consider itinerancy and carrier super- 
fluidity. Within this lowest lobe we may take hardcore a 
and b bosons. This will also be convenient for the numer- 
ical simulations in section [V] Whilst this does not affect 
physics within the lobe, the zero hopping upper boundary 
is modified by the restriction on the 6-atom population. 
In this case we only need retain the states |0,0), |1,1) 
and the admixtures of |0, 1) and |1, 0). The zero hopping 
diagram shown in Fig. [2] is replaced by Fig. [4] The lower 
boundary remains unchanged because the same eigen- 
states are involved, and the upper boundary becomes a 
mirror reflection of the lower one. This is quite natu- 
ral since the a and b operators now appear on an equal 
footing, modulo the effects of the band splitting. More 
formally, this may be traced to the invariance of the hard- 
core Hamiltonian ([5]) (up to a constant term) under the 
particle-hole transformation a — >• a', b — >• b\ together 
with fix — > e a + £b + U a b — Mi, and the interchange of the 
a and b operators and the hopping parameters J a and 
Jb', since hardcore bosons obey onsite anticommutation 
relations, n — > 1 — n, under particle-hole transformation. 
For e a = — 1, eb = 1 and U a b = this involves fXi —> —jJL\. 
The particle-hole transformation also accounts for the 
change in the vacuum, 1 0, 0) — > |1, 1). 
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FIG. 4: (Color online). Variational zero hopping diagram for 
hardcore a and b atoms. The symmetry about /ii = reflects 
the combined particle-hole and species interchange symmetry 
of the Hamiltonian (JS]) in the hardcore limit; see text. 



To incorporate itinerancy and superfluidity we aug- 
ment the variational analysis for two component bosons 
in an optical lattice HU with a coherent state for light: 



|V> = |7> ® [cos0i(cosxiaj + smxib\) 



+ sin 9i(cos i]i + simjiblaj) |0) 



(34) 



where (7) is the coherent state introduced previously, and 
Xj Vi 7 are to be determined. The corresponding or- 
der parameters are given by (a) = 4 sin 20 cos (x — rj), 
(b) = Isin2#sin(x + rj) and (-0) = 7. The first term in 
brackets in equation (fMl) describes the Mott state, and 
the second superfluidity. For = this coincides with 
the variational approach for localized excitons coupled to 
light [33j, and as we will discuss in section HVBl repro- 
duces the previous results for J a = 0. More generally, 
(|34[) takes real hopping into account, involving site va- 
cancies and interspecies double occupation. It provides a 
useful starting point to identify the boundaries between 
the Mott and superfluid regions. We consider spatially 
uniform phases with energy density £ = (V\H\V) /N: 

£ = (e + — e_ cos 2%) cos 2 9 + (2e + + U a b) sin 2 77 sin 2 9 

- I [J a cos 2 (x -f]) + J b sin 2 (x + rf)\ sin 2 29 

+ ujj 2 + 97 cos 2 9 sin 2%, 

where z is the coordination and e± = (e& ± e a )/2. Min- 
imizing on 7 gives 7 = (ip}/y/N = — g cos 2 9 sin 2x/2w. 
This may be eliminated from (|35|) to yield 

£ =(e+ - e_ cos2x) cos 2 6> + (2e + + f7 ah ) sin 2 ?/sin 2 9 

- - [J a cos 2 (x -rj) + Jb sin 2 (x + ??)] sin 2 29 



g 2 cos 4 sin 2 2\ 




4a3 



(36) 



FIG. 5: (Color online). Variational phase diagram with J a = 
Jb = J and e a = —1, ej, = 1, ui — g — 1, U a b = 0. The phases 
are (i) a Mott insulator (MI, dark blue), (ii) a superradiant 
Mott state supporting a condensate of photo-excitations (SR- 
MI, light blue), (iii) a superradiant superfluid (SR-SF, light 
red), and (iv) an a-type superfluid (a-SF, dark red). As shown 
by the solid line, the tetracritical line ultimately bifurcates 
into two bicritical points connected by a first order transition. 
For > 0, the a-type superfluid is replaced by a b-type owing 
to the particle-hole and species interchange symmetry of the 
hardcore Hamiltonian ([5]) when J a = Jf, see text. 



The symmetries of (|36[) determine the domain of mini- 
mization. It is invariant under 9 — > 9 + n and 9 — > —9, 
and may be minimized over 9 G [0,7r/2]. Likewise, it is 
invariant under \ ~~ ^ X + 71 " an d f] — > T)+ir. Since the hop- 
ping contribution favors x an d f] taking the same sign, we 
may also minimize x> V S [0, tt/2] . To begin with we set 
Ja = Jb = J and the expression may be further reduced 
using J a cos 2 {x~ r l) + Jb sin 2 (x + v) = J(l + sin2xsin2?7). 
Minimizing over the restricted domain yields the phase 
diagram in Fig. [5] For the chosen parameters, we have 
up to four distinct phases in the interval (i\ < 0; (i) a 
Mott state with (a) = (b) = (tp^ip) = 0, (ii) a superra- 
diant Mott state with (a) = (b) = and (ip^ip) ^ 0, 
(iii) a single component superfluid with (a) 7^ and 
(b) = (ijj^ijj) = 0, and (iv) a superradiant superfluid 
(a) ^ 0, (b) 7^ 0, {iftij)) ^ 0. Indeed, the Hamiltonian 
displays aU(l)xU(l) symmetry and these may be broken 
independently. The phase diagram reflects this pattern of 
symmetry breaking. In particular, the superradiant Mott 
state corresponds to an unbroken U (1) in the matter sec- 
tor (corresponding to a pinned density and phase fluctu- 
ations) but a broken U(l) (or phase coherent condensate) 
for photo-excitations. As we shall discuss below, the ex- 
pectation value of the bilinear, (b^a) ^ 0, corresponds to 
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the onset of coherence in the Dicke model. This novel 
phase may thus be regarded as a form of supersolid [62j 
in which photo-excitations condense on the background 
of a Mott insulator. Although the lattice precludes spon- 
taneous translational symmetry breaking (at least with 
this periodicity) the excitations may be thought of as 
mobile defects in an otherwise ordered background. 

As in the zero hopping case, we may extend Fig. [5] into 
the region /ii > by exploiting symmetries of ([5]). This 
is reflected in the variational energy by using sin 2 n = 
(1 — cos2rj)/2 to combine the e+ contributions: 



„ I. . 1 . 



2 [Uab - (e a 



£6 



1 — ^2) cos 2\ cos 2 9 
U a b — 2/xi) cos 2rj\ sin 2 9 



zJ 



( 1 + sin 2\ sin 2ry) sin 2 26> 



<7 2 cos 4 1 



! 2x 



4(w - /Lt2) 



(37) 

Since J a = Jt = J this is invariant (up to a constant) 
under /ii — s- e a + ej, + U a b — Mi an d *7 — > 7r /2 — which 
interchanges the superfluid order parameters (a) and (6) . 
For the parameters used in Fig. [5l the a-type superfluid 
observed for fix < is replaced with a 6- type for fix > 0. 

As may be seen from Fig. [5j the locus of the tetracrit- 
ical line may be determined from the intersection of the 
superradiance transition with the Mott insulator to a- 
type superfluid phase boundary. This yields e a + z J = 
as discussed in AppcndixO This may also be seen from a 
Landau type expansion of £. In view of the non- linear re- 
lationship, it is convenient to expand in the angles 9, Xi V 
as opposed to the order parameters (a) , (b) , (ip) : 



£ = e a - {e a + zJ)9 2 



UJLJq 



(38) 



The quadratic "mass" terms vanish when e a + zJ = and 
the superradiance condition g = ^JujCjq is met. In some 
simple cases £ may be expanded directly in the order 
parameters. For example, throughout the entire Mott 
lobe where (a) = (b) = we may set 9 = 0. Returning 
to equation (|35[) and minimizing over x yields tan 2\ = 
-2<77/w . This gives £ = e + + ujj 2 - y/ufi/A + g 2 */ 2 , 
where 7 = (ip) / \N. This agrees with our variational 
zero hopping result (fTT|) in the lowest lobe with n = 1. 

In Fig. [S] we present a cross section of the phase dia- 
gram ([5]) for (ii = —0.6. We indicate continuous transi- 
tions by single lines and first order transitions by dotted 
lines. In general, the transition from the superradiant 
Mott state to the non-superradiant a-type superfluid in- 
volves a discontinuous jump in the photon density. This 
results in a first order transition as indicated by the dis- 
continuity in the first derivative of the energy and the 
order parameters as shown in Fig. [7] This first order line 
connects two bicritical points as indicated in Fig. [6] The 
length of this first order segment changes with hi, and 
emanates from the tetracritical line as shown in Fig. [S] 
These overall features are also exhibited for more general 
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FIG. 6: (Color online). Cross section of the phase diagram 
J5j for /ii = —0.6. The first order transition from the super- 
radiant Mott state to the a-type superfluid is indicated by a 
dashed line. The locus of this transition is determined analyt- 
ically from (|C8p . The remaining transitions are continuous. 



hopping parameters as shown in Fig. [HI where a-type and 
6-type superfluids emerge for large hopping asymmetry. 



B. Coherence in Superradiant Mott Phase 

An interesting aspect of the superradiant Mott phase is 
that polariton condensation coexists with the Mott char- 
acter. Within both Mott phases 9 = and our variational 
wavefunction (|34[) becomes 

IV) = |7) ® II ( cos X^l + smxibf) |0). (39) 

i 

To describe photo-excitations above the filled Mott state, 
we introduce a change of vacuum \Q) = rii a l|0> so that 

|V) = |7>®n cos ^n (l + tanxWo*) ( 40 ) 



Since we are dealing with hardcore bosons this may be 
exponentiated, and for homogenous parameters 

|V) = I7) g> e Nhlcosx e tanx ^ b * a t\n). (41) 

This is already reminiscent of a coherent state for the 
bilinears, although one needs to be careful since we are 
dealing with hardcore bosons. Instead, we may examine 
condensation properties directly by computing expecta- 
tion values using (|39[) . Using the Schwingcr boson repre- 
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FIG. 7: Evolution of the mean field order parameters (a), 
(b), 7 = (ijj) /\/N across the transitions depicted in Fig. [6] for 
fi2 = 0.45. (I) First order transition from the superradiant 
Mott state to the a-type superfluid, (II) continuous transition 
from the a-type superfluid to the superradiant superfluid. In- 
set: Evolution of the variational energy (relative to the zero 
hopping Mott phase) across the transitions in Fig. [6] for fixed 
values of ^2- The solid line at /12 = 0.45 has a discontinuity 
in the first derivative and indicates a first order transition in 
passing from the superradiant Mott insulator to the a-type 
superfluid. The remaining transitions are continuous. 
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FIG. 8: (Color online). Variational phase diagram showing a 
slice through Fig. [4] with fii = —0.25 and extended in to the 
(J a , Jb) hopping plane. We use the same key as in Fig. [5] and 
denote the 6-type superfluid (b-SF) in orange. Single compo- 
nent superfluids are observed for large hopping asymmetry. 



sentation 



for the lowest lobe 



NUMERICAL SIMULATIONS 



we may calculate the bosonic bilinear (J + ): 



nt-nf), 



(42) 



(V\J+\V) = 2<V|6k|V) = jVsinxcosx = E s in2 X 



(43) 



W+|V) _ E^k 



N 



N 



(44) 



This agrees with the result ([55)1 obtained from the 
Holstein-Prirnakoff approach to the Dicke model [59| . 
We see that in the two-band Bose-Hubbard problem we 
have condensation involving particle-hole pairs above the 
Mott background with (b'a) 7^ 0. This corresponds to 
condensation of the Holstein-Prirnakoff boson (c) 7^ 
in the dual formulation as evidenced by equation (|33|) . 
Crucially, condensation of (wa) is not accompanied by 
condensation of the individual bosonic carriers a and b. 
This follows immediately from (f39| where (a) = (b) = 0. 
Since matter-light coherence in the Dicke model trans- 
lates into polariton and not carrier condensation, Mott 
behavior and polariton condensation may coexist. 



We now analyze the Hamiltonian (JTJ by exact diag- 
onalization, imposing a maximum number of photons, 
(ip'ip) < M^, in addition to the hardcore a, b con- 
straints. We consider a one-dimensional system with 
N lattice sites, with the basis of tensor product states 

\<t>) { a a) ® \<t>t b) ® \<t>)^ ] « - - 1 ON 



where v a b = 1, • • • , 2 and 



0,1,..., M, 



V.- 



Here 



<y a ) 

L a,l > 



L a,N/ 



where n y "™' e {0, 1}, and the photon states are given 



At 6 = our variational analysis yields 7 = — g sin 2x/2cl>: ^ v 



lip 



where r&* ] e {0, 1, 



.,M^}. The 



total Hilbert space has dimension D = 2 N 2 N (M^ + 1). 
In the largest case considered, with TV = 8 sites and 
AI,,p = 64 photons, this corresponds to a matrix dimen- 
sion D « AE6. Increasing further has only minor 
influence. 

Applying periodic boundary conditions in real space, 
we diagonalize the sparse matrix representation of H to 
obtain the groundstate |$o) and its energy Eq. We com- 
pute the atom and photon densities, and the density fluc- 
tuations a a = yj (n„) — (h a ) 2 . To obtain the superfluid 
fraction, of the a atoms, we impose a phase twist 
<C 7T by means of a Peierls factor a\a^ i-> a\a^e~ 1@ / N , 
and calculate the change in the ground state energy [63| : 



N E, 



En 



J a. (^a) 



e 2 



(45) 



For the single-component Bose-Hubbard model, this 
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quantity is zero deep in the Mott insulator, and ap- 
proaches unity far in the superfiuid. Note that in our 
case, the hardcore constraint always provides an effec- 
tive interaction even for large hopping, so that f s < 1. A 
total superfiuid density can be obtained by imposing the 
phases on both species a,b and calculating 



N 



02 



e 



(46) 



We supplement these superfiuid diagnostics with the zero 
momentum occupations n a (k = 0) = iV _1 ^2r>„(ata q ) . 
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A. Small Hopping and Dicke Superradiance 

To verify the existence of a superradiance transition 
within the Mott phase we begin our numerical treatment 
in the limit of small non-zero hopping. Figure [9] shows 
the Mott lobe with density n = 1 at z J = 0.1, for two dif- 
ferent values of M$ . For sufficiently large M$ the results 
are in good agreement with the variational zero-hopping 
phase diagram shown in Fig. [4j and the overall features 
depend only weakly on the number of sites. In order to 
examine the onset of superradiance, we track the evolu- 
tion of the polariton and photon densities in Fig. [TUJ At 
zero hopping, equation (|24|l yields the critical coupling 
~g\ = uuq = (u> — /i2)(eb — (a — ^2) ■ For the chosen pa- 
rameters at fixed 9=1, the transition occurs at a critical 
chemical potential, = (3 — \/5)/2 « 0.382. This onset 
is well reproduced at small finite hopping, as shown in 
Fig- QUI For N = 1 (where the Dicke model reduces to the 
Jaynes-Cummings model) we see quantized steps with 
fixed integer polariton density, N%/N. These form the 
basis of the Mott lobes observed in Jayncs-Cummings- 
Hubbard models [9j, [23|. Even in this extreme limit, the 
densities closely track the thermodynamic Dicke model 
results. As N increases, the step sizes are reduced by 
a factor of 1/N and we approach the variational ther- 
modynamic results. This behavior in our itinerant bo- 
son model closely mirrors direct finite N simulations of 
the Dicke model; see Fig. [16] A notable difference be- 
tween the Dicke and Jaynes-Cummings-Hubbard mod- 
els [H, [23[ (both with N sites) is reflected in their zero 
hopping cigenstatcs. In the latter, eigenstates are tensor 
products of superpositions of two states which differ in 
photon number by one, whereas in the Dicke model the 
global photon mode leads to a coherent photon state for 
large N. In the Dicke model the higher Mott lobes are 
eliminated in favor of a continuous photon onset, whereas 
in the Jaynes-Cummings-Hubbard models [t| [23| they 
remain in tact. 



B. General Phase Diagram 

Turning to the construction of the overall phase dia- 
gram, we first consider a fixed value //i = and present 



FIG. 9: Exact diagonalization results for the phase diagram at 
small hopping, zj — 0.1, for hardcore a,b atoms and e a = — 1, 
£(, = 1, u = g = 1, U a b = 0. The dashed (solid) lines corre- 
sponds to iV = 8 and = 16 (M^ = 64) photons, and have 
been obtained from spline fits to the contours at n = 0.99 and 
n = 1.01. The vertical line indicates the location of the super- 
radiance transition as determined from (ip 1 tp)/N > 0.01. The 
results are in good qualitative agreement with the variational 
zero hopping results in Fig. [4] showing the continuity of this 
behavior into the Mott phase. 
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FIG. 10: Exact diagonalization results with N = 1,8 sites, 
J a — Jb = J and zj = 0.1. We set e a = —1, e b = 1, 
uj = g = 1, hi = Uab = 0. The profiles show the onset 
of superradiance within the Mott phase, and evolve with in- 
creasing system size towards those for the thermodynamic 
limit of the Dicke model. Evolution of (a) the photon density 
{ip^tp) /N, (b) the polariton density, N 2 /N. For N = 1 the 
Dicke model reduces to the Jaynes-Cummings model and ex- 
hibits quantized integer steps in the polariton density |23(|. 
More generally these are quantized in units of 1/N. The value 
of M^, is of minor importance in the range of \ii considered as 
it sets the upper limit on the number of excitations available. 
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results in the (zJ,^) plane; sec Figures [TT1 and [T2| So 
far we have set U a b = in order to expose the essential 
details. Following our original work [361 ] . we now take 
U a b = 1 to help illustrate the generality of the overall 
results. The interspecies interaction also helps stabilize 
the Mott region up to larger values of the hopping. 

In Fig. QT] we show density plots of atom and photon 
densities, and density fluctuations. In Fig. Q2] we plot 
supcrfluid fractions. From these data, we can identify 
the four different phases present in Fig. [SJ The Mott- 
supcrfluid transition is visible from the deviation of the 
total atom density from unity in Fig. [TT](a). The su- 
perradiance transition corresponds to the onset of pho- 
ton density in Fig. [TT](b). Panels (c) and (d) demon- 
strate the existence of the single-component supcrfluid 
in the region where n ^ 1 and n a ^= 0, while n& = 0. 
This is also visible in the individual superfluid fractions 
shown in Fig. [12] We also see the onset of a non-trivial 
population imbalance, nb — n a , or existence of 6-bosons, 
accompanying the supcrradiance transition. Specifically, 
n a = 1, rib = in the left of the white region of panel (b), 
but n a < 1, rib > above; see panels (c), (d) and (g). The 
nature of the two Mott phases (normal and superradiant) 
is illustrated by panels (e)-(f). We observe an onset of 
local fluctuations a a with increasing hopping. Since the 
Mott state is characterized by n = n a + n& = 1, the indi- 
vidual atom density fluctuations a a in the superradiant 
Mott phase are much larger than a. In the normal Mott 
phase, rib = so that the Mott insulator consists purely 
of a-atoms and hence a a = <Jb = 0. Finally, the Mott- 
supcrfluid transition is also apparent in our results for 
the superfluid fractions as shown in Fig. [T21 

Exact diagonalization of small clusters yields an ap- 
proximation to the critical value J c for the Mott- 
supcrfluid transition by monitoring the onset of f s , as in 
Fig. [T2](a), or the local fluctuations a of the total atom 
density, as in Fig. [TT (e) . Both a and f s are also nonzero 
below J c for finite N 63| . While f s is expected to scale to 
zero for J < J c as N — > oo, the local density fluctuations 
remain nonzero in the Mott phase due to virtual hopping 
processes. This is illustrated in Fig. Q2] The system size 
N has a noticeable impact on the onset of correlations re- 
lated to superfluidity. Since a finite size scaling analysis 
requires much larger system sizes, we provide approxi- 
mate phase boundaries obtained from the contour lines 
o{zJ,\i2) = 0.3 (alternatively f s (zJ, [12) = const.) to in- 
dicate the onset of superfluidity, and (ip'ip)/N = 0.01 
to indicate superradiance. The numerical constants for 
the contours have been chosen in order to obtain phase 
boundaries that match those suggested by the data in 
Figs. ITT1 and [T2l see Fig.[l3j Figure[T3](c) shows the zero- 
momentum occupation nik = 0) = n a (k = 0)+n&(fc = 0) 
for the atoms. The latter is expected to diverge as a func- 
tion of system size in the SF phase, and such behavior 
can indeed be seen above z J ps 1.1. 

Repeating the above procedure for different \i\ values 
we build up a picture of the overall phase diagram as 
shown in Fig. 1141 This may be compared to the varia- 
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FIG. 11: Exact diagonalization results for N = 8 sites and 
M,/, = 16 photons. We set /ii = 0, J a = Jb = J and e a — — 1, 
e& = 1, u) = g = Uab = 1. The panels show density plots of 
(a) total atom density, (b) photon density, (c) a-atom density, 
(d) b-atom density, (e) fluctuations of the total atom density, 
(f) fluctuations of the a-atom density, (g) fluctuations of the 
b-atom density. In particular (a) reveals the Mott-superfluid 
transition and (b) the superradiance transition. 
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FIG. 12: Exact diagonalization results for N = 8 sites and 
— 16 photons. We set = 0, J a = Jb — J and e a = —1, 
eb = 1, w = g = Uab = 1. The panels show density plots 
of (a) total atomic superfluid fraction, (b) a-atom superfluid 
fraction, (c) b-atom superfluid fraction. The Mott-superfluid 
transition is evident from the onset in panel (a), and the single 
component a-type superfluid in panels (b) and (c). 



tional analysis shown in Fig. [5j Wc find the same phases 
as in the analytical approach and the evolution of the 
phase boundaries is in good agreement. For the choice 
of nx = these phases meet in a tetracritical point. As 
found analytically this extends into a tetracritical line 
which ultimately bifurcates into two bicritical points. 
The agreement between the numerical simulations and 
mean field theory is remarkable given the enhanced role 
of fluctuations in low dimensions. This mirrors the suc- 
cess of mean field theory in other bosonic systems and 
may be assisted by the long range cavity photons. 



VI. RELATION TO OTHER PROBLEMS 

A feature not addressed by the present mean field the- 
ory, but captured in Fig. 1141 is the dispersion of the su- 
perradiance transition with J; in the Mott phase, 6 = 0, 
and J drops out of the variational energy (pj|) . One way 
to understand this is to recast the matter contribution as 
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+ sinx4a i )|ft), (47) 



is the filled Mott state in the ab- 



where |f2) = Yii 
sence of excitations. In the spatially uniform case 

|V M ) =Me x ^ b ^\n), 



0.8 



0.6 



0.4 



0.2 



(a) Superfluid fraction 




N=4 






N=6 






N=8 


■-'/ 






JzJ c -1.1 




(c) Zero-momentum 
occupation 




(48) 



O 3 
II 

~c 2 
1 



FIG. 13: Exact diagonalization results with J a = Jb = J 
and e a = — 1, et, = 1, ui = g = U a b = 1. We set fii = 0, 
fi2 — 0.4, M^j = 2N, and show the system size dependence 
of (a) superfluid fraction, (b) density fluctuations, (c) zero 
momentum occupation. Within the Mott phase for J < J c , 
the superfluid fraction decreases with increasing system size. 
This indicates a Mott-superfluid transition at approximately 
zj c m 1.1. This is compatible with the fluctuation onset cri- 
terion, a = 0.3, that we use to determine the overall phase 
diagram in Fig. 1141 and the divergence of n(k = 0) with sys- 
tem size in panel (c). 



where A = tan% and J\f = (1 + A 2 )~ w/2 . This only ac- 
commodates local particle- hole pairs above the filled Mott 
background, and may thus be regarded as the BEC limit 
of the BEC-BCS crossover problem. By analogy with the 
fermionic BCS approach to cxciton insulators [HI, Ir35| . 
and the crossover problem in 40 K [66j], one expects that 
excitons may lower their energy by spreading out in real 
space and pairing in momentum space. In this regard, the 
"BCS" pairing phenomenon in Bose gases has a long his- 
tory, with the Valatin-Butler wavefunction (67| playing 
the role of the BCS state. This has been developed in a 
series of early works motivated by liquid 4 He (68l - [73j and 
biexciton formation (74T - t76j . This later emerged in stud- 
ies of High- T^i superconductivity [77[ and two-component 
Bose gases [78|, [79[ . The interesting possibility of bound 
states of three or four particles has also been explored 
[80l . Hl[ . Central to these studies is the requirement to 
stabilize attractive Bose gases against collapse, with the 
aid of short-range repulsion or internal structure. When 
this condition is met, and at sufficiently low densities 
where exciton overlap is negligible, it has been argued 
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FIG. 14: (Color online). Overall phase diagram obtained by 
exact diagonalization with N = 8 sites and = 16 photons. 
We set J a — Jt = J, £a = —1, e& = 1, oj = g = 1, and U a b = 1, 
and use the same key as in Fig. [5] We extract the Mott- 
superfluid boundaries from the onset of density fluctuations, 
a > 0.3, and the superradiance transition from the onset of 
photons, {ipiip)/N > 0.01; see text. The simulations show 
the distinct phases and the bifurcation of the tetracritical line 
suggested by mean field theory. The bifurcation is indicated 
by the regions where the two boundaries overlap along a line 
rather than a single tetracritical point. 

that such pairing states may exist (73|. In the present 
context, stabilization of the paired state and the absence 
of carrier condensation is brought about by the interplay 
of Mott physics and photoexcitation. It would be in- 
teresting to explore this problem in more detail, and we 
leave these refinements to future studies. 

As noted in our previous work, the connection to the 
BEC-BCS crossover for bosons is reinforced by the Fes- 
hbach resonance problem studied in the continuum (821 - 
l86j and on the lattice [87l - l89j . Performing a particle-hole 
transformation, the matter-light coupling reads ip mbi- 
Aside from the global nature of the photon, this converts 
a and b into a "molecule" ip. At the outset there are 
eight possible phases corresponding to separate conden- 
sation of (a) , (6) , (ip) . At the mean field level, only five of 
these may survive; condensation of two variables provides 
an effective field (as dictated by the coupling) which in- 
duces condensation of the other. In Figures 151 and fl4l the 
band asymmetry, e a < ef,, in conjunction with the cho- 
sen parameters, reduces this to four or less. Nonetheless, 
the additional b-type superfluid is stabilized for larger 




a b 

FIG. 15: (Color online.) The classical light limit of equation 
([T]) (where ip is a replaced by a c- number), may be simulated 
in optical superlattices where the coupling g represents 
tunnelling between different wells, a and b. 



values of fi± owing to the particle-hole and species in- 
terchange symmetry involving fii — > e a + €(, + U a h —JM ■ 
In contrast to the single species mean field theory [82| - 
[Hj], this two species case supports an atomic superfluid, 
since condensation of one carrier no longer induces an 
effective field. Moreover, condensation may leave a U(l) 
symmetry intact, which allows the coexistence of Mott 
and phase coherent behavior. 

In deriving (|35[) and the phase diagram, we are primar- 
ily concerned with the matter-light coupling. As such we 
incorporate U a h as in Ref. l56l . Within this variational ap- 
proach, this interaction has no impact on the Mott states 
with 9 = 0, as evident from the energy density (|35p . 
Nonetheless, the presence of the matter-light coupling 
stabilizes the non-trivial phases in Fig. [5] and provides 
good agreement with the numerical simulations. How- 
ever, as noted by Soyler et al [9(| and previous works [78h - 
Hl|, analogous pairing phases may be supported in the 
two-component Bose-Hubbard model, without matter- 
light coupling, through a more sophisticated treatment of 
U a b itself. Indeed, onsitc repulsive interactions, U a bn a nb, 
favor a particle of one species and a hole of the other on 
the same site. Treating this pairing in a BCS approach, 
one may replace n°n- by |Ai| 2 + (A^fa^ + h.c), where 
Aj = {a\bi), is to be determined self-consistently. This 
field acts as a local "photon" , and a similar mean field 
phenomenology may ensue. Such pairing also occurs in 
fermionic models [37|. Although our discussion has fo- 
cused on a single global photon, the symmetry analysis 
is more general. This is supported by studies of the two- 
band Bose-Hubbard model for equal fillings and com- 
mensurate densities (9lj . 

In closing we note that the classical light limit of equa- 
tion |T]) (where ip is replaced by a c-number) may be 
simulated in optical superlattices [92| where giaib\ rep- 
resents tunnelling between different wells; see Fig. [131 In 
the case of hardcore a-atoms and softcore 6-atoms this 
provides an analogue of the Jaynes-Cummings-Hubbard 
models considered in Refs. [23|. This geometry may 
also be useful in realizing other "matter-bath" systems. 
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VII. CONCLUSIONS 

We have investigated the two-band Bose-Hubbard 
model coupled to a single mode of a cavity light field. 
We combine analytical and numerical techniques and find 
good agreement between the approaches. The model dis- 
plays a novel phase in which "polaritons" condense on 
the platform of a bosonic Mott insulator. We extend 
our previous work }3r| in several directions including an 
investigation of the overall phase diagram, and the na- 
ture of polariton condensation. In particular, we use 
the framework of the Dicke model to discuss how po- 
lariton condensation emerges in the absence of carrier 
condensation. In terms of numerical results, we have 
presented superfluid fractions, atom density fluctuations 
and zero-momentum occupation, and analyzed the phase 
diagram under variation of In addition, we have 
addressed the effects of finite cluster size and the pho- 
ton cutoff. This helps illustrate connections to work on 
Jaynes-Cummings-Hubbard models and coupled cavity 
arrays. This topic has broad connections to other prob- 
lems of current interest including atom-molecule mix- 
tures and the BEC-BCS crossover in bosonic systems. 

There are many avenues for further research includ- 
ing non-equilibrium aspects and collective excitations 
[93I |94| . It would be interesting to develop numerical 
techniques to explore finite temperature polariton con- 
densation and the onset of phase coherence in the Mott 
phase. It would also be worthwhile to examine the phase 
diagram with softcore bosons with finite U aa and Ubb- 
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Appendix A: Dicke Models 



written in the so-called rotating wave approximation in 
which terms of the form Sfip^ are excluded. In the ab- 



sence of these terms, the model is integrable [42J, |4J, |58[ 
and exhibits a superradiance quantum phase transition 
when g = -Jujljq [4ll444j | ; for a discussion of the model in- 
cluding counter- rotating terms see Ref. 54 1. Depending 



on the application these models are often recast in differ- 
ent representations and we gather a few results below. 



1. Jaynes— Cummings 

The simplest case is the Dicke model with N = 1, or 
the Jaynes-Cummings model, describing a single two- 
level system coupled to radiation. It plays a central role 
in the variational analysis presented in section fill Al and 
recent problems in matter-light systems, e.g. [9|, [23J. A 
convenient representation is the single site Hamiltonian 



H 



-ji a + e b nb + g(b^a + alb), 



(A2) 



where a is a hardcore boson (restricted to occupancy zero 
or one) and b may take arbitrary occupancy. We may 
construct exact eigenstates of total integer occupation, 



n > 1, as superpositions 



In) = all, n — 1) 



(A3) 



where the first and second entries of the states are the oc- 
cupations of the a and b particles respectively. In matrix 
form this yields the eigenvalue problem 



e a + (n - l)e b g^/ri 
g^fri ne b 



(A4) 



where we recall that b\n) = \/n\n — 1) and b^n) = 
\Jn-\- \\n + 1) for bosons. The eigenvalues are given by 



E 



ne b — wq/2 ± ^/wq/4 + g 2 n, 



(A5) 



Throughout this manuscript we make use of the Dicke 



model 41 



and its reductions. In the literature this ap- 
pears in various guises and with different names so we 
present a brief guide. The Dicke model describes N two- 
level systems or "spins" coupled to radiation 



JV 



9_ 

N 



N 



H = cj J2 S i+ "V'V + T^Sr + Sfi>), (Al) 



where Si is a spin-1/2 operator and ip is a canonical bo- 
son. Here ujq is the level-splitting between the two-level 
systems, uj is the frequency of the cavity mode, and g 
is the strength of the matter-light coupling. In quan- 
tum optics the Dicke model is also known as the Tavis- 
Cummings model [§5j]. In the special case of a single 
spin, N = 1, the Dicke model is often referred to as the 
Jaynes-Cummings model [47j. The Hamiltonian (|A1[) is 



where loq = e b — e a . This agrees with the spectrum of 
the Jaynes-Cummings model given by Greentree et al 0] , 
where we identify e a , e b , and g with their parameters e, uj, 
and (3. The coefficients in the corresponding normalized 
eigenvectors \n}± may be written 



a± 



At 



-oj /2±x{n) 



(A6) 



where x( n ) = + 9 2n - We note that there is a mi- 

nor typing error in the eigenstates given in the methods 
section of Ref. 0, where the basis states |l,n — 1) and 
|0,n) are erroneously reversed. Equivalcntly, one may 
reverse the sign of A = ljq in their expression. 
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FIG. 16: Evolution of the polariton density, N2/N, in the 
Dicke model for N = 1, 2 sites and comparison to the TV — > 00 
thermodynamic limit. For JV = 1 the Dicke model reduces 
to the Jaynes-Cummings model and the quantized polariton 
steps correspond to the Mott lobes discussed in Ref. In 
the Dicke model these steps are quantized in units of I/TV and 
as JV -> 00 wc approach the variational results. Inset: Evolu- 
tion of the photon density. The magnetization (or population 
imbalance) may be obtained by subtraction. All figures show 
the onset of superradiance at [12 = (3 — y/E)/2 ~ 0.382, for 
u — ¥ Hi = 1 — /i2, Wo — >• uio = 2 — (1,2, and g = 1. 



Dicke Model 



The generic Dicke model (|A1|) is integrable for arbi- 
trary N [ill, HH, (Hll , but is most conveniently analyzed 
in the thermodynamic limit, N — > 00 [59| . However, in 
numerical simulations we must deal with the effects of 
finite N and truncations of the photon Hilbert space. In 
Fig. [16] we show the evolution of the polariton density, 
N 2 /N, where N 2 = iftij) + X^ 5 ? + 1/2), and the photon 
density as a function of N. This is obtained by solution 
of the finite dimensional matrix problem. Both plots ex- 
hibit discrete jumps which track the thermodynamic re- 
sults. For N = 1 the Dicke model reduces to the Jaynes- 
Cummings model and the quantized polariton steps cor- 
respond to the Mott lobes discussed in Ref. 



Appendix B: Zero Hopping Phase Diagram with 

3 = 



An instructive way to think about the topology of the 
zero hopping phase diagram in Fig. [2l is in the absence 
of the matter-light coupling. In Fig. [T7] we plot the loci 
e a = 0, If, = 0, uj = 0, corresponding to population 
transitions in the Hamiltonian ((6|) for g = 0. When g is 
switched on, the horizontal boundaries bend downwards 
and continuously evolve into those shown in Fig. [5J 



FIG. 17: Zero hopping phase diagram for the Hamiltonian ([5]) 
in the absence of matter-light interaction, with e a = — 1, ej, = 
1, uj = 1, g = 0. The loci e a = 0, q, = 0, uj — 0, indicated by 
the lower, upper and vertical lines respectively, delineate the 
regions of population onset for a, b, ip in the grand canonical 
ensemble. When g is switched on the horizontal boundaries 
bend downwards and evolve into those shown in Fig. [2] 



Appendix C: Variational Phase Boundaries 

The variational energy for the finite hopping problem 
allows some analytic progress with the phase boundaries, 
and highlights connections to the Dicke model. In the 
absence of competition from other phases, the transition 
between the non-superradiant insulator (8 = x — 7 — 0) 
and the a-type superfluid (9 7^ 0, \ = V = 7 = 0) f° r 
example occurs when e Q + zJ = 0. This may be seen by 
explicit computation of the energies of each phase. The 
energy density in the generic Mott phase is 



£e=o = (e+ - ? - cos2x) 



f sin 2 2 X 



(CI) 



Minimizing on \ yields either x — 0, corresponding to an 
ordinary Mott insulator, with £mi = ~ £- 1 or 



cos 2x 



UlUJQ 



(C2) 



corresponding to the superradiant Mott insulator, 
this non-trivial stationary point 



At 



d 2 £ 
dx 2 



\9 ^0) 
Cog 2 



(C3) 



This corresponds to a minimum (with Co > 0) provided 
g > \/ujujq. This coincides with the superradiance tran- 
sition in the Dicke model [4H - I44I ] . The energy density in 
the superradiant Mott phase is 



(g 2 -wwo) 2 



-SRMI 



'MI 



(C4) 



Similarly, in the a-type superfluid 



zJ 



£ aSF = (?+ - e_ ) cos 2 6 - ^ sin 2 20. (C5) 
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Minimizing on 9 yields either 9 = or cos 29 = — £ ^ 
The latter yields 



zJ 



aSF 



(£+ - g- - zJ) s 
4zJ 



(C6) 



The condition £mi = £aSF yields the MI-aSF boundary 
e + — 6- + 2 J = or equivalently e a + zJ = 0. This meets 
the superradiancc onset g 2 = wwq at g 2 = + z J) or 

Vi = e b - ^--J^ + Z J. (C7) 

2 W - /Z 2 

In a similar fashion, the transition between the supcrra- 
diant Mott state and the a-type superfluid occurs when 

^SRMI = £aSF Or 



zJ(g -QQ ) =gu>(e a +zJ) 



(C8) 



Appendix D: Absence of Higher Lobes in Coupling 
to Quantum Light 

In Fig. [2] we see that only the Mott lobe with atomic 
density, n = 1, is supported due to the coupling to the 
quantum light held. To prove that only this lobe survives 
we compare the energies of the higher states to those de- 
picted in Fig. Without loss of generality we may focus 
on the stable region defined by q, — g 2 /4cj > 0. Within 



this domain E~(-f vllI ) and E~(j = 0) are both smoothly 
increasing functions of n. Furthermore, whilst the super- 
radiance condition (fl4")l holds, E~(0) > £?~(7 var ). Wc 
immediately conclude that whenever the n = 1 superra- 
diance condition is met (i? 4 > uj 2 uj 2 ), the ground state 
energy is cither E^{-y vax ) or 0, the vacuum energy. How- 
ever, outside of this region we must examine whether the 
lowest energy supcrradiant state is lower in energy than 
either E± (0) or the vacuum. 

To address this let us consider superradiant states 
with n > 1. These are candidates for the ground state 
within the region, g 4 < ui 2 Qq < n 2 g 4 , where the up- 
per bound follows from the superradiance condition and 
the lower bound precludes the region in which -Ej~(7 V ar) 
is supported. If ujq > 0, taking roots of the inequal- 
ity yields <2>o > uiuj 2 /ng 2 and ljq > g 2 /Q. In con- 
junction with the condition for stability, we see that 
£v7(7var) > e a = E^(0). That is to say, the ground 
state never has n > 1. 

If &o < a stronger statement can be made. Invoking 
the superradiance condition, ujq < —uiu 2 lng 2 , it follows 
that E~(j var ) > for all n. Hence the vacuum must 
be the ground state when Qq < 0. (The latter implies 
~ ?a > £fc > 0, where the last inequality has 
followed from the stability requirement). We see that a 
sufficient condition for termination of the Mott lobe is 
fj,2 > ljq. However, since ^2 < w in general, such regimes 
are only stable for ujq < u; see Fig. O 
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